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g : ABSTRACT 

o : 

CN ■ 

' We calculate the reionization history in Cold Dark Matter (CDM) models. The 

, ^ ■ epoch of the end of reionization and the Thomson scattering optical depth to the cosmic 

microwave background depend on the power spectrum amplitude on small scales and 
on the ionizing photon emissivity per unit mass in collapsed halos. We calibrate the 
emissivity to reproduce the measured ionizing background intensity at z = 4. Models in 
Tlj- I which all CDM halos have either a constant emissivity or a constant energy emitted per 

' Hubble time, per unit mass, predict that reionization ends near z ~ 6 and the optical 

■ depth is in the range 0.05 < Te < 0.09, consistent with WMAP results at the la to 2(t 

^ ■ level. If the optical depth is as high as 0.17 (as suggested by WMAP), halos of velocity 

. dispersion 3 — 30kms~^ at z > 15 must have ionizing emissivities per unit mass 

, larger by a factor > 50 compared to the more massive halos that produce the ionizing 

emissivity at z = 4. This factor increases to 100 if the CDM power spectrum amplitude 
O \ is required to agree with the Croft et al. (2002) measurement from the Lya forest. If 

' Tg > 0.17 were confirmed, a higher ionizing emissivity at z > 15 compared to z = 4 

might arise from an enhanced star formation rate or quasar abundance per unit mass 
^ ■ and an increased escape fraction for ionizing photons; the end of reionization could have 

^ . been delayed to z ~ 6 because of the suppression of gas accretion and star formation in 

^ . low-mass halos as the medium was reionized. 

Subject headings: cosmology: theory - diffuse radiation - intergalactic medium - galax- 
ies: formation 



1. INTRODUCTION 

The epoch of reionization of the universe started with the emission of the first ionizing photons 
into the intergalactic medium (hereafter, IGM), and ended when all the low-density regions of the 
universe were ionized. During this period of time, the ionizing photons reaching the IGM had to 
be sufficient to ionize every atom in the universe and to balance any recombinations in the IGM (in 
addition, many more ionizing photons may have been emitted in regions of dense, self-shielded gas, 
which were locally absorbed and did not contribute to ionizing the IGM). At the end of reionization, 



- 2 - 



the mean free path of ionizing photons increased to scales much larger than the typical separation 
of the collapsed halos hosting the sources. Regions of dense, self-shielded gas (observed as Lyman 
limit systems) shrank in size as the ionizing background increased in intensity, due to the growing 
mean free path (Miralda-Escude, Haehnelt, &: Rees 2000). If the fraction of baryons in regions 
of intermediate density which were ionized as the mean free path increased was small, the rise in 
the background intensity could be relatively fast, as seen in the numerical simulations of Gnedin 
(2000). 

Precisely this rapid change of the ionizing background intensity has been inferred by Fan et 
al. (2002) to occur at 2; ~ 6, by analyzing the change of the mean transmitted flux of the Lya 
forest in the highest redshift quasars known. The transmitted flux apparently drops abruptly to 
very low levels near z = 6. This observation most likely implies that the end of reionization 
occurred, in fact, at 2; ~ 6, although as a larger number of quasars at 2; > 6 are discovered, more 
detailed investigations of this issue will need to be done to understand precisely how the mean free 
path grew. In particular, the number of small gaps of transmitted flux, such as the one seen in a 
z = 6.37 quasar by White et al. (2003), provides information on the size of the Hii regions before 
they overlapped (Miralda-Escude 1998; Barkana 2002). 

The results of the Wilkinson Microwave Anisotropy Probe ( WMAP; Bennett et al. 2003 and 
references therein) have strengthened the observational support of the Cold Dark Matter (CDM) 
scenario in a flat model that contains a vacuum energy component (in addition to the ordinary 
baryonic matter and dark matter). The model parameters are being measured with increasing 
accuracy via combined observations of the cosmic microwave background (CMB), galaxy clustering, 
weak lensing, and the Lya forest. The theory then makes an increasingly reliable prediction for 
the abundance of collapsed dark matter halos as a function of mass and redshift. It is in these 
halos that the process of radiative gas dissipation and star formation can take place, leading to 
the emission of ionizing radiation. However, the uncertainties in the efficiency of star formation, 
the stellar initial mass function (IMF), the emission due to gas accretion into black holes (be they 
in X-ray binaries or active galactic nuclei), and the fraction of the ionizing photons that escape 
into the IGM from the dense regions in which they are emitted, prevent a clear prediction for the 
history of reionization. 

Even if the time at which reionization ends is determined to be at z 2± 6 by the Gunn-Peterson 
trough (Gunn & Peterson 1965) observations, there can be a long time interval in which the universe 
was partially ionized, with cosmological Hii regions around the sources filling only some fraction 
of the volume. A very good observational probe for the characteristic time at which the bulk of 
the baryons in the universe were ionized is the electron scattering optical depth to the CMB, Tg. 
This has been measured for the first time by WMAP (Kogut et al. 2003; Spergel et al. 2003) using 
the large-scale polarization-temperature anisotropy correlation. At present, the result still has a 
large error, both observational and model-dependent: Kogut et al. find a model-independent value 
of Te = 0.16 lb 0.04 from the polarization-temperature correlation of the WMAP data, but Spergel 
et al. find Tg = 0.17 it 0.06 from a fit of a CDM model with a running spectral index to a combined 
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data set of CMB, galaxy clustering and Lya forest observations, and = 0.12 ± 0.06 when the 
fit is forced to the usual CDM model with a primordial power-law spectrum. A value of Tg = 0.17 
would imply that most baryons were reionized as early as z ~ 20. 

The announcement of the measurement of Tg by WMAP has led to a large number of papers 
discussing the possibility of this early reionization in the CDM model (Whyithe &: Locb 2003; 
Haiman k. Holder 2003; Holder et al. 2003; Ciardi, Ferrara, & White 2003; Somerville & Livio 
2003; Sokasian et al. 2003; Cen 2003; Chiu, Fan, & Ostriker 2003). These authors have generally 
found that, if one is willing to assume a high efficiency in the production of ionizing photons per 
unit mass in the low-mass halos collapsing at early times, then early reionization and a high optical 
depth is possible in CDM models. Naturally, the lower the power spectrum amplitude on small 
scales, the higher the efficiency that is required, and some models can be ruled out based on their 
small power on small scales, such as the warm dark matter model (Spergel et al. 2003). But to a 
large extent, the ability of the CDM model to accommodate a wide range for the epoch at which the 
bulk of the IGM was reionized simply reflects our ignorance of the efficiency in producing ionizing 
photons that escape. 

However, a reasonable question we can ask is how the required efficiency compares to any 
observationally determined values of the rate at which ionizing photons have been emitted over 
the history of the universe, taking into account the predicted abundances of halos of different 
masses that form at each redshift. In this paper, we adopt a very simple model of the ionizing 
emissivity from CDM halos, but which still captures the essential physical ingredient for the process 
of reionization: as the medium is ionized, the gas is heated and the Jeans mass is raised, implying 
that some low-mass halos where gas can collapse when the medium is neutral are unable to continue 
to accrete gas after the medium is ionized. We divide halos into two populations, the "high-mass" 
ones (which will be called Population A) that arc able to accrete gas and form new stars irrespective 
of the state of the intergalactic medium, and the "low-mass" ones (Population B) that can form 
stars and emit ionizing radiation only when the medium is neutral. Low-mass halos collapse first 
in CDM, but as they ionize the medium around them the rate of star formation can be suppressed, 
until the high-mass halos start to collapse (Couchman & Rees 1986; Shapiro, Giroux, k. Babul 
1994). A similar separation of halo classes was made by Haiman & Holder (2003); our Population 
A corresponds to their Type II and Type la halos (Haiman & Holder separated these halos into 
the ones that can or cannot cool by atomic cooling after they have accreted the gas; we will not 
distinguish between these two classes in our more simple model), and our Population B corresponds 
to their Type lb population. 

Our calibration of the efficiency to produce ionizing photons is based on observations of the Lya 
forest during the post-reionization era. The mean transmitted flux of the Lya forest, combined with 
CDM simulations of Lya forest spectra and measurements of the baryon density of the universe, 
allows us to infer the intensity of the ionizing background; and the abundance of Lyman limit 
systems tells us the mean free path of ionizing photons. The ratio of the background intensity to 
the mean free path then yields the volume-averaged emissivity. This emissivity was inferred to be. 
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at most, seven ionizing photons per baryon and per Hubble time at 2: = 4 by Miralda-Escude (2003). 

Remaining uncertainties on the value of this emissivity are related mainly to the temperature of the 
IGM and the CDM power spectrum amplitude. At z > 4, the abundance of Lyman limit systems 
becomes highly uncertain, so the emissivity cannot be inferred. The key question for reionization 
is, of course, how this emissivity changes with increasing redshift. If the comoving emissivity stays 
constant at z > 4, the optical depth Tg is not larger than 0.09 (Miralda-Escude 2003). In CDM 
models, the fraction of mass in halos that can form stars decreases with redshift ai z > 4, so to 
obtain Tg > 0.09 (which requires an increase of the emissivity with redshift), the low-mass halos 
(or Population B) formed at z ^ 4 must have a much higher efficiency to emit ionizing radiation 
than the halos of higher mass at ^ ~ 4. 

The question we are addressing in this paper is: how much higher does this efficiency have to 
be in order to produce a certain value of the optical depth Tg, and what constraints are imposed by 
requiring reionization to end at 2; = 6? We describe our model for the efficiency to emit ionizing 
radiation in §2. Results are presented in §3, and a general discussion is given in §4. 

2. MODELS 

2.1. CDM Power Spectrum 

We use three CDM models to predict the halo abundances. All the models assume flat space 
with a cosmological constant, with Q^nO = 0.27, Q^q = 0.044, h = 0.71. The three models vary in 
the normalization and slope of the power spectrum: Model 1 has as = 0.84, Ug = 0.93; Model 2 has 
as = 0.9, Ug = 0.96; and Model 3 has ag = 0.9, Ug = 1 (here, as is the linear rms mass fluctuation 
at present on spheres of radius 8h~^ Mpc, and Ug is the primordial scalar spectral index). The 
power spectrum is computed with the fltting formula of Eisenstein &: Hu (1999). Model 1 has the 
smallest amplitude of fluctuations on the small scales on which the first halos that emit ionizing 
radiation collapse (0.03 to 1 Mpc), and Model 3 has the largest amplitude. Model 3, with the 
largest power on small scales, has halos forming at the highest redshift, and therefore predicts an 
earlier start of reionization and higher optical depth Tg than the other models for a fixed ionizing 
emissivity per unit mass. 

Model 1 has the same parameters as the best fit model found by WMAP (Bennett et al. 2003; 
Spergel et al. 2003) to the combined data set of CMB observations, the 2dF Galaxy Redshift Survey 
(Colless et al. 2001), and Lya forest power spectrum data (Croft et al. 2002). However, we do not 
include the running spectral index introduced in Spergel et al. (inclusion of the running spectral 
index would reduce the power on small scales). Model 2 has higher values of as and Ug by 1.5 
and 1 standard deviations, respectively, from the best fit values of Spergel et al. Model 2 also has 
a slightly larger power than the best fit CDM model without a running spectral index found by 
Spergel et al. (which has parameters cjg = 0.8 and Ug = 0.96). Finally, Model 3 is very close to the 
best fit of Spergel et al. to CMB data alone, and is the same model that was adopted in various 
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recent studies on the reionization history (Ciardi et al. 2003; Somerville & Livio 2003; Sokasian et 
al. 2003; Cen 2003). 

Measurement of the Lja forest has allowed a determination of the power spectrum amplitude 

on the smallest scales. We therefore compare the amplitude of our models with the determination of 
the Lya forest power spectrum by Croft et al. (2002), who find a normalization at redshift z = 2.72 
of 

A'(M^^^(M = 0.741°:?°, (1) 

where H{z)~^kp = 0.03(km/s)~^. The normalization of our Models (1, 2, 3) are A'^{kp) = 
(0.85,1.06,1.27), respectively. Hence, Models 2 and 3 have a power spectrum amplitude higher 
than the Croft et al. (2002) measurement by l.Qa and 2.6a", respectively. The high amplitude of 
Model 3 would seem to be ruled out by the result of Croft et al.; note, however, that the quoted 
error on this measurement does not include some potentially important systematic effects, as dis- 
cussed in Croft et al., and the amplitude can be higher if the mean Lya transmitted flux is higher 
than the value used by Croft et al. (see Seljak, McDonald, & Makarov 2003). 



2.2. Halo Abundances 

We use the Sheth & Tormen (1999) prescription for computing the abundances of halos as 
a function of mass and redshift. Halos are divided into two populations: those which do not 
emit any ionizing photons after the medium around them has been reionized (Population B), and 
those which are massive enough to continue forming stars and producing ionizing photons after 
reionization (Population A). In this paper we place the division between the two populations at a 
velocity dispersion ctq = 35kms~^, based on the suppression of gas infall and cooling in smaller 
halos caused by photoionization (Efstathiou 1992; Thoul &: Weinberg 1996; Navarro & Steinmetz 
1997). The lower limit to the velocity dispersion of the Population B of low-mass halos is placed 
at (^min = 3.68 km s~^, because molecular hydrogen cooling below a temperature of ~ 2000 K is 
ineffective (Yoshida et al. 2003). 

The fraction of baryons, Fa,b, which occupy each of the halo populations is shown in Figure 1 
for Models 1 to 3. We also show in Figure 2 the difference between using the Sheth-Tormen and the 
Press-Schechter (Press & Schechter 1974) mass functions, for Model 1. The Sheth-Tormen mass 
function predicts a slightly larger number of objects at the highest redshifts. 
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2.3. Emissivity 

For Population A, i.e., halos with a > ctq = 35kms^^, the ionizing photon emissivity per 
baryon is assumed to be of the form 

..(.)=^.W^(i|^)", (2) 

where Fa{z) is the fraction of baryons in Population A halos at rcdshift z, and £4 is the total 
emissivity at z = A, which is assumed to be equal to 7 ionizing photons per baryon and per Hubble 
time at 2; = 4, the upper limit obtained in Miralda-Escude (2003). We will use two values of a in 
this paper: for a = the energy emitted per unit physical time and per unit mass remains constant, 
and for a = 1.5 the total energy emitted per Hubble time at redshift z remains constant. Note that 
gas cooling is inefficient in very massive halos when the virialized gas temperature is > 10^ K, but 
at z = 4 the abundance of these halos, which arise from > 4a fluctuations, are negligibly small, so 
our assumption of a constant emissivity per unit mass for Population A is reasonable. 

For Population B, we allow an emissivity per unit mass different from the one in Population 
A halos, parameterizing it as 

where 7 is a scaling factor used to adjust the variation in efficiency between the low-mass and 
high-mass halos. This results in a smooth variation of the emissivity with a. The Population B 
emissivity is then 

eeiz) = / fB{(J,z)e{a)da , (4) 
where fs is the fraction of mass in halos in the interval da around a, and we set (7 rain — 

3.68kms"\ 

corresponding to 2000 K for the mean particle mass of the primordial atomic gas. 



2.4. Analytical Framework 

The equation for the ionization of the intergalactic medium is (Madau, Haardt, &; Rees 1999; 
Miralda-Escude 2003) 

= eA + esil -y) - Ry, (6) 

where y is the volume fraction of the IGM that is ionized at any given time t, and R is the 
mean number of recombinations per baryon in the ionized regions per unit time. This equation 
makes the approximation that the photon mean free path is very short, so emitted photons are 
instantaneously used to ionize atoms. It also hides all the complicating physics of reionization in the 
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averaged recombination rate R, which depends on the clumping factor of the photoionized gas. Note 
that Population B halos emit ionizing radiation only when they are located in the atomic medium, 
so the emissivity is multiplied by 1 — y. We assume a clumping factor of unity throughout this 
paper, computing the recombination coefficient at T = 10^ K, which gives R = 1.16 X 10-1'^ s-i at 
z = 4. The solution of equation (6) is 

y{t)= Jle^p[-fl[R{t")+eB{t")]dt"} 

X [eB{t') + eA{t')] dt' , (7) 

where U is an arbitrarily chosen initial time (before any appreciable halo formation). 



2.5. Optical Depth 

The Thomson scattering optical depth, Tg, is 

^ nefJeC r Fi{l + zfy{z) dz 
''"^'^ Ho Jo [OAo + f^™o(l + ^)3]i/2 ' ^ ^ 

where ^l\o = 1 — ^mO, n-e is the comoving electron number density, (Tg is the electron cross section, 
and Fi is the fraction of the baryons in the ionized regions that are actually in the ionized IGM 
(with the rest being in self-shielded atomic or molecular clouds, or in stars). We assume = 0.9 
(as in Miralda-Escude 2003) for z < 6, but Fj = 1 at z > 6; lowering Fj at z > 6 would demand 
starting reionization at higher redshift to obtain a fixed optical depth, increasing the required 
emissivity from the first star-forming halos (note that equation (5) needs to be modified if Fi < 1 
during reionization). For our value of Q,t,h^ = 0.022, the electron number density is ng = 2.7 x 
10~^(1 — fyy) cm~^, where Y = 0.24 is the helium abundance by mass, and /y = (0.5, 0.75) when 
the helium is (doubly, singly) ionized. We assume that helium becomes doubly ionized at z = 3.5 
(e.g., Heap et al. 2000), and that it is singly ionized at higher redshift wherever hydrogen is ionized. 
This results in an optical depth up to z = 6 of Te{z = 6) = 0.032. 



3. RESULTS 

We now present the results on the fraction of the volume that is ionized as a function of 
redshift, y{z), with the emissivity that arises from our two populations of halos given by equations 
(2) and (3). To summarize the emissivity model, the Population A (high-mass halos) emissivity 
per unit mass is fixed at z = A and is allowed to vary as (1 -|- z)", while the Population B (low-mass 
halo) emissivity per unit mass can also be increased with the parameter 7. 

Figure 3 shows y{z) for each of the three CDM models and for both a = and a = 1.5. The 
case 7 = 0, where the emissivity per unit mass from Population A and B is the same, is shown 
as solid lines. The cases 7 = 1, 10 and 100 are shown as dotted, dashed, and dash-dot lines. 
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respectively. The models with high 7 values have an early reionization as the first Population B 

halos collapse; then, y{z) grows slowly when it approaches unity because of the decrease in the 
emission from Population B as most of the IGM is ionized, and finally y reaches unity when the 
Population A halos start forming and emitting radiation in the ionized regions of the universe. 

The optical depth to electron scattering for all the models shown is listed in Table 1. 

The cmissivity as a function of redshift is shown in Figure 4 for the same models. At redshift 
just above 4, the emissivity declines with redshift because the abundance of Population A halos 
(which are the only ones contributing when the medium is entirely ionized) drops rapidly. The 
emissivity increases again at redshifts higher than the end of reionization, when Population B halos 
contribute. Models with higher 7 can obviously have a higher emissivity at very high redshift, but 
their emissivity declines at intermediate redshifts when y is close to unity and the emission from 
Population B halos is suppressed, as assumed in equation (3). The models with higher 7 result in 
higher optical depths because the IGM is reionized earlier. 

4. DISCUSSION 

For a fixed CDM halo population, the most simple model for the emissivity is to say that the 
rate of photon emission per unit mass is constant for all halos in which cooling can take place, i.e., 
with a > CFmin- This model, using the emissivity per unit mass calibrated with the observational 
determination at z = 4 and a clumping factor of one, predicts a redshift at which reionization ends 
of 5 to 7, in very good agreement with the observation of the appearance of the Gunn-Peterson 
trough at z = 6. The emissivity history in this model shows a decline with redshift from 2 = 4 
to 6, then a slight increase up to z ~ 10 as more Population B halos contribute, and then a rapid 
decline beyond z ~ 10 to 15 as the number of all halos decreases. 

Obviously, the reionization history should also depend on the clumping factor affecting the 
recombination rate. We have assumed a clumping factor of unity, and a higher clumping factor 
would delay the end of reionization. A moderate clumping factor of 2 or 3 would not greatly affect 
the result, because the recombination time at the mean baryonic density is equal to the Hubble 
time at z = 6 and the increase in y{z) is relatively fast for the 7 = models (solid lines in Fig. 
3; this is simply due to the rapid increase in the mass fraction in collapsed halos shown in Fig. 
1), so the number of recombinations is not very large. For example, a clumping factor of 3 for 
Model 1 with a = 1.5 (panel b) delays the end of reionization from z = 6 to z = 5. As discussed 
by Miralda-Escude et al. (2000), and Miralda-Escude (2003), the clumping factor is unlikely to be 
large during reionization. Therefore, we have the remarkable result that the simplest model of a 
constant emissivity per unit mass in CDM halos, calibrated with the observed emissivity at z = 4, 
roughly predicts the right epoch for the end of reionization. 

This same model predicts an optical depth to Thomson scattering, Tg, of 0.05 to 0.09, depending 
on a and on the allowed variations on the power spectrum parameters determining the amplitude 
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of fluctuations on the small scales of the flrst halos. Although this is smaller than the value found 
by the WMAP mission in a CDM model fit according to Spcrgel et al. (2003), Tg = 0.17 ± 0.06, 
it is not ruled out with the present errorbar: the upper value of Tg = 0.09 of our 7 = models 
is consistent with this measurement at the 1.3a level. The likelihood function of this range of 
Te shown in Figure 8 of Spergel et al. (2003) is also seen to be not very small. Moreover, if the 
running spectral index in the CDM model is not allowed, the best fit then yields Tg = 0.12 it 0.06 
(see Table 7 of Spergel et al.). The smaller errorbar obtained by the analysis in Kogut et al. (2003) 
is puzzling, because the addition of observational data and the use of a parameterized model instead 
of a model-independent analysis should not result in an increased error. 

If a high value of Tg is confirmed, then the ionizing emission efficiency in the earliest low-mass 
halos to form in CDM is required to be much higher than in the sources present in the more massive 
halos at z = 4. The required factor by which the efficiency must increase can be inferred from Table 
1. The increase in emissivity per unit mass is ~ (1 + 27)(1 + 2)° (see eq. 3), for a typical halo 
contributing to the emissivity at high redshift with a ~ Skms"^. For the case a = 0, the factor 
1 + 27 gives the increase in the rate of emission per unit mass, while for the case a = 1.5, the factor 
1 + 27 gives the increase in the total energy per unit mass emitted in ionizing photons, taking into 
account the shorter Hubble time at higher redshift. Interpolating from Table 1, we find that to 
reach Tg = 0.17, the total ionizing energy per unit mass emitted by Population B halos at 2; > 15 
needs to be larger by a factor 1 + 27 ~ 10 compared to the Population A halos at 2; = 4, and the 
rate of emission needs to be higher by a factor ^ 100. For Model 1 (more consistent with the Croft 
et al. 2002 measurement of the Lya forest power spectrum), the energy emitted per unit mass needs 
to increase by a factor ~ 25, and the emission rate per unit mass needs to increase by ~ 300. 

Are these very high eflEiciencies for producing ionizing photons in the earliest low-mass halos 

physically realistic? The emissivity of Population A halos at z = 4 is about 60 photons per baryon 
per Hubble time, as obtained by dividing the mean emissivity £4 of 7 photons per baryon per Hubble 
time (see §2.3) by the fraction of baryons in Population A halos at z = 4 (see Fig. 1). From the 
emission ratios mentioned above required to produce Tg = 0.17, Population B halos at 2; ~ 20 would 
need to emit ~ 1500 and 600 photons per Hubble time for Models 1, and 3, respectively. Since 
in these models with high optical depth the universe is rcionized at z ~ 20 over a time interval 
corresponding to Az ~ 5 (see Fig. 3), or one quarter of a Hubble time, the total emission during 
this time interval coming from the first halos needs to be about 400 and 150 photons per baryon 
in Population B halos, for Models 1 and 3. 

To compare this to the number of photons that can be produced by star formation, we use 
the models of Leitherer et al. (1999; see their Fig. 78) to find that for a Salpeter IMF with solar 
metallicity in the stellar mass range IMq < M < IOOM0, 6000 ionizing photons per baryon are 
emitted. Most of the ionizing photons are produced by stars with M > 20Mq, and in the case of a 
Salpeter IMF with stars forming only in the range 20AfQ < M < WOMq, 25000 photons per baryon 
can be emitted. Note that this number is equal to the ratio of the energy obtained by hydrogen 
fusion (7 MeV) to the average energy of an ionizing photon (20 eV) times the fraction of stellar mass 
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that is fused to helium (~ 30 %) times the fraction of energy emitted as ionizing photons (~ 25%). 

Very massive stars with zero metalHcity can further increase this emission to 10^ photons per 
baryon (Tumhnson & Shull 2000; Bromm, Kudritzki, &: Loeb 2001; Schaerer 2002), because their 
convective cores during the main-sequence include almost all the mass of the star (zero metallicity 
stars also have higher temperature than metal-rich ones during the main-sequence, which increases 
the fraction of energy they emit in ionizing radiation, but also increases the average energy of their 
ionizing photons). Thus, if is the fraction of the gas in a halo that can be converted to stars with 
M > 2OM0, and /esc is the fraction of ionizing photons that escape to the IGM, then an optical 
depth Tg = 0.17 requires /*/esc > (10~^'^, 10~^'^) for Models (1,3) for zero-metallicity stars, and 
f*fesc > (10^^-^, 10^^-^) for Models (1,3) for metal-enriched stars. 

It is clear from these values that an optical depth Tg = 0.17 can be produced with physically 
realistic models. Even Tg = 0.3 is possible in Model 3 with 7 = 100 and a = 1.5, which implies 
an emission of 3000 ionizing photons per baryon (or ^^.^f^sc = lO"-*^-^ for zero-metallicity stars) from 
the first star-forming halos over the redshift range 25 < z < 32. However, if a high optical depth 
is confirmed by future measurements, an explanation for the dramatic change in ionizing emission 
from the low-mass halos that formed the first stars to the more massive halos at lower redshift will 
need to be found. 

A difficulty in making star formation very efficient in low-mass halos is that their shallow 
potential wells makes them more vulnerable to lose their g result of stellar winds, pho- 

toionization and supernovae. The escape fraction might be higher than in more massive halos at 
lower redshift, but only after gas has been ionized and expelled from the vicinity of stars, in which 
case the gas can easily escape from the halo. The typical halo where stars form at z 2± 20 with 
velocity dispersion ~ 5kms~'^ and total mass of 10^ M0 (e.g., Yoshida et al. 2003) has a gas mass 
of 2 X 10^ Mq within a virialized radius ~ 200 parsecs, and cooling by molecular hydrogen leads 
to the formation of a dense central gas clump containing only ~ 1% of the gas mass. The rest of 
the gas must take more than a halo dynamical time (~ lO'''^ years) to reach the center and cool 
before it has any possibility to form stars, by which time supernovae would have taken place. The 
hydrodynamic simulations of these first star-forming halos (Abel, Bryan, & Norman 2002; Bromm, 
Coppi, & Larson 2002) suggest that the first stars may all be very massive. If the ~ 2OOOM0 of 
the core of cooled gas can efficiently form one or a few massive stars, and most of the emitted 
ionizing photons escape, one may barely be able to achieve the fraction /*/esc ~ 10~^'^ required for 
Te = 0.17. After the formation of these stars, the supernovae explosions would expel the remaining 
gas from the halo (Bromm, Yoshida, & Hernquist 2003). 

In halos of higher mass (with a of 10 to 30 kms^^), a larger number of supernovae would be 
required to expel the gas. Nevertheless, the effects of photoionization and supernova explosions 
should still self-regulate the rate of massive star formation, and it is not clear why the earliest 
galaxies at high redshift could have higher values of f^fesc than the more massive galaxies formed 
at lower redshift. While the mean gas density in galaxies should increase with redshift, increasing 
gas cooling rates, the shallower potential wells of low-mass halos allow for gas dispersal after a 



- 11 - 



smaller fraction /* of gas has turned to massive stars. 

Another possibility is that black holes formed by the first massive stars are able to accrete gas. 
Gravitational accretion can produce radiation with an efficiency ~ 20 times greater than stellar 
fusion (and stars generally fuse only a fraction of their mass). In this case, requiring a fraction as 
low as ~ 10^^'^ of the halo gas to accrete on a central black hole could produce the '-^ 300 photons 
per baryon in halos required for Tg = 0.17, depending on the fraction of the accretion energy emitted 
as ionizing photons (we have assumed 20% here). However, whether one speculates about a high 
efficiency of massive star formation or of gas accretion into black holes, the question that remains 
unanswered is why such high efficiency to produce ionizing radiation would occur only in the first 
halos at high redshift, and not in similar or more massive halos at later times. 

If it were nevertheless possible to produce these high emissivities when the first stars formed, a 
second possible problem of an early reionization might be that the high emissivity leads to a rapid 
end of reionization, in contradiction with the presence of the Gunn-Peterson trough at z = 6. Our 
results suggest a way to solve this particular problem: as seen in Figure 3, once the fraction of the 
IGM that is ionized becomes close to unity, the emissivity can be decreased because star formation 
in Population B halos is suppressed. This can regulate the rate of emission to a level that just 
balances recombinations, and it is only when the more massive Population A halos reach a critical 
abundance that the emissivity can rise again, completing reionization and causing the fast increase 
in the mean free path of ionizing photons that clears the Gunn-Pctcrson trough. Thus, in all our 
models the end of reionization occurs between redshift 5 and 8, even though the models with high 
7 can ionize half of all the volume as early as 2; ~ 20. In the models with the earliest reionization 
in Figure 3, the redshift for the end of reionization could easily be lowered to fit the observed value 
of 6 by moderately increasing the clumping factor after y is close to unity or increasing the value 
of (To that separates our Populations A and B. 

Note that models with high 7 have an emissivity history that is double peaked: there is a first 
maximum due to Population B halos, then the emissivity declines as most of the IGM is ionized, 
and a second maximum is then caused by Population A halos. The possibility of a doublc-pcakcd 
history of star formation due to the effect of reionization raising the Jeans mass was proposed by 
Couchman Sz Rees (1986) and Shapiro et al. (1994). However, this does not necessarily imply a 
"double reionization" , as other authors have argued (Cen 2003; Whyithe &; Loeb 2003). The fraction 
of ionized volume y{z) rises monotonically with time, and it is only "stalled" by the suppression 
of the star formation due to ionization, in a self-regulated process. A double reionization in which 
the universe recombines again after having been ionized would require a long delay of the feedback 
ionization effects on the star formation rate. If this delay were present, the suppression of star 
formation might start to operate only after stars have already emitted more than enough photons 
to ionize all the universe; later, when the IGM recombines again, star formation might not resume 
at the same level for some time, allowing for partial recombination to occur. This feedback delay 
might be caused by a reservoir of gas in Population B halos that can keep forming stars after the 
surrounding low-density medium is ionized, and by the finite time required for new recombined and 
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cooled gas to accrete in the Population B halos after the IGM recombines again; but it is not clear 
if this would be sufficient to cause the global instability in the star formation rate that would be 
required for a double reionization. 

Finally, comparing our results with other recent work, we note that other authors (e.g., Whyithe 
& Locb 2003; Haiman &: Holder 2003; Ciardi ct al. 2003; Somerville & Livio 2003; Sokasian ct al. 
2003; Cen 2003) have generally emphasized that the central value of the H/7lf^P result Tg 0.17 can 
be obtained in CDM models with a physically allowed emissivity of ionizing photons. Although 
we agree with this conclusion, we note that a drastic change in the emissivity per unit mass of 
star-forming halos from z > 15 to z = 4 is required. For example, Ciardi et al. (2003) find, using 
our Model 3, that Tg = 0.16 can be reached in models where the dominant emission comes from 
halos corresponding roughly to our Population A, but the emissivity they assume would vastly 
overpredict the intensity of the ionizing background at z = 4. At the same time, models with no 
such change in the emissivity yield a lower optical depth that is still within the plausible error 
range of the WMAP result. Our approach is closest to that of Chiu et al. (2003), who use an 
ionizing emission efficiency chosen to yield reasonable models for the ionizing background, and our 
conclusions are generally in agreement with theirs. 



5. CONCLUSION 

The first detection of the Thompson optical depth to the CMB by the WMAP mission intro- 
duces a new era in the study of reionization. The observations of the Gunn-Petcrson absorption 
in the highest redshift quasars, which are sensitive to the time when reionization ended and the 
mean free path of ionizing photons increased, can be confronted with the total ionized column to 
the surface of last scattering, which is sensitive to the epoch when about half of the baryons in 
the universe were ionized. The simplest assumption we can make about the emissivity is that it is 
roughly constant in all halos in which gas can cool. When we calibrate the emissivity to the upper 
limit allowed by observations at ^ = 4, we find that the end of reionization is close to z = 6, as 
suggested by observations, but that the optical depth is smaller than 0.09. To obtain an optical 
depth as high as Tg = 0.17 (the central value of the WMAP result), an increase of the rate of ion- 
izing emission per unit mass in collapsed halos with redshift, of a factor ~ 100 to 300 (depending 
on the CDM power spectrum model) from z = 4 to z ~ 20, is required. 

Further progress on our understanding of the reionization epoch will come, among other things, 
from an increased accuracy in the measurement of Tg, as well as more robust measurements of the 
power spectrum on small scales from Lya forest observations. The discovery of quasars at z > 6 
and the search for any small fraction of Lya transmitted fiux in gaps along their Gunn-Peterson 
trough (White et al. 2003) may also allow us to test if the fraction of the ionized medium was close 
to unity for a wide redshift range above z = 6 (as in our models with high 7 in Fig. 3), or if it 
increased rapidly over a short redshift range. If an optical depth as high as the central value of 
the WMAP measurement is confirmed, then either a physical explanation for the high efficiency 
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of the ionizing emissivity of the first objects will need to be found, or some more fundamental 
modification of the CDM model, such as non-Gaussian fluctuations on small scales (e.g., Avelino 
k. Liddle 2003), may be implied. 

This work was supported in part by NSF grant NSF-0098515. 

REFERENCES 

Abel, T., Bryan, G. L., & Norman, M. L. 2002, Sci, 295, 93 

Avelino, P. P., & Liddle, A. R. 2003, MNRAS, submitted (astro-ph/0305357) 

Barkana, R. 2002, NewA, 7, 85 

Bennett, C. L., et al. 2003, ApJ, accepted (astro-ph/0302207) 
Bromm, V., Coppi, P. S., & Larson, R. B. 2002, ApJ, 564, 23 
Bromm, V., Kudritzki, R. P., & Loeb, A. 2001, ApJ, 552, 464 

Bromm, V., Yoshida, N., & Hernquist, L. 2003, ApJ, submitted (astro-ph/0305333) 
Gen, R. 2003, ApJ, submitted (astro-ph/0303236) 

Ghiu, W. A., Fan, X., & Ostriker, J. P. 2003, ApJ, submitted (astro-ph/0304234) 

Ciardi, B., Fcrrara, A., & White, S. D. M. 2003, MNRAS, submitted (astro-ph/0302451) 

Golless, M., et al. 2001, MNRAS, 328, 1039 

Gouchman, H. M. P., & Rees, M. J. 1986, MNRAS, 221, 53 

Groft, R. A. G., Weinberg, D. H., Bolter, M., Buries, S., Hernquist, L., Katz, N., Kirkman, D., & 
Tytler, D. 2002, ApJ, 581, 20 

Efstathiou, G. 1992, MNRAS, 256, 43 

Eisenstein, D. J., & Hu, W. 1999, ApJ, 511, 5 

Fan, X., Narayanan, V. K., Strauss, M. A., White, R. L., Becker, R. H., Pentericci, L., k, Rix, 
H.-W. 2002, AJ, 123, 1247 

Gnedin, N. 2000, ApJ, 535, 530 

Gunn, J. E., & Peterson, B. A. 1965, ApJ, 142, 1633 

Haiman, Z., & Holder, G. P. 2003, ApJ, submitted (astro-ph/0302403) 



-14- 



Heap, S. R., Williger, G. M., Smette, A., Hubeny, I., Sahu, M., Jenkins, E. B., Tripp, T. M., 
Winkler, J. N. 2000, ApJ, 534, 69 

Holder, G., Haiman, Z., Kaplinghat, M., & Knox, L. 2003, ApJ, submitted (astro-pli/0302404) 

Kogut, A., et al. 2003, ApJ, accepted (astro-ph/0302213) 

Leitherer, C, et al. 1999, ApJS, 123, 3 

Madau, P., Haardt, F., & Rees, M. J. 1999, ApJ, 514, 648 

Miralda-Escude, J. 1998, ApJ, 501, 15 

Miralda-Escude, J. 2003, ApJ, submitted (astro-pli/0211071) 
Miralda-Escude, J., Haehnelt, M., & Rees, M. J., ApJ, 530, 1 
Navarro, J., & Stcinmctz, M. 1997, ApJ, 478, 13 
Press, W. H., & Schechter, P. 1974, ApJ, 187, 425 
Schaerer, D. 2002, A& A, 382, 28 

Seljak, U., McDonald, P., & Makarov, A. 2003, MNRAS, accepted (astro-ph/0302571) 
Shapiro, P. R., Giroux, M. L., & Babul, A. 1994, ApJ, 427. 25 
Sheth, R. K., & Tormen, G. 1999, MNRAS, 308, 119 

Sokasian, A., Abel, T., Hernquist, L., & Springel, V. 2003, ApJ, submitted (astro-pli/0303098) 

Somerville, R. S., & Livio, M. 2003, ApJ, submitted (astro-pli/0303017) 

Spergel, D. N., et al. 2003, ApJ, accepted (astro-pli/0302209) 

Thoul, A. A., & Weinberg, D. H. 1996, ApJ, 465, 608 

Tumlinson, J., k Shull, J. M. 2000, ApJ, 528, L65 

White, R. L, Becker, R. H., Fan, X., & Strauss, M. A. 2003, AJ, 126, 1 

Whyithe, S., & Loeb, A. 2003, ApJ, submitted (astro-ph/0302297) 

Yoshida, N., Abel, T., Hernquist, L., & Sugiyama, N. 2003, ApJ, in press (astro-ph/0301645) 



This preprint was prepared with the AAS lAT^^X macros v5.0. 



-15- 




Redshift 

Fig. 1. — Fraction of mass collapsed into halos as a function of rcdshift for two ranges in halo 
velocity dispersion, a. Population A includes the halos with a > ao, while Population B represents 
the halos with amin < a < ao- For each of the Model 1 (solid), Model 2 (dashed), and Model 3 
(dotted) cosmological parameter sets, amin = 3.68 km and ao = 35kms~^. 



-16- 




Fig. 2. — Comparison of the amount of collapsed mass in halo Populations A and B as a function 
of redshift derived with the Sheth-Tormen (solid) and Press-Schechter (dashed) mass functions for 
Model 1. 
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20 30 10 
Redshift 

Fig. 3. — Reionization history of the universe for six cases: (a) Model 1, a = 0; (h) Model 1, 
a = 1.5; (c) Model 2, a = 0; (d) Model 2, a = 1.5; (e) Model 3, a = 0; (f) Model 3, a = 1.5. 
Within each panel, the lines indicate the reionization history for 7 = (solid), 1 (dotted), 10 
(dashed), and 100 (dot-dash). 
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Redshift 

Fig. 4. — Log emissivity (ionizing photons per baryon per Hubble Time at z = 4) versus redshift. 
Panels and lines as in Figure 3 



Table 1. Model Optical Depths 



Model 1 Model 2 Model 3 



7 Te{a = 0) Te{a = 1.5) Te(a = 0) Te{a = 1.5) Te{a = 0) Te(a = 1.5) 






0.047 


0.071 


0.050 


0.080 


0.053 


0.088 
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0.060 


0.092 


0.064 
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0.068 


0.117 


10 


0.098 


0.146 


0.108 


0.171 


0.117 


0.194 


100 
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